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While Josephson-junction-like structures intrinsic to the layered cuprate high temperature super- 
conductors offer an attractive stage for exploiting possible applications to new quantum technologies, 
the low energy quasiparticle excitations characteristically present in these d-wave superconductors 
may easily destruct the coherence required. Here we demonstrate for the first time the feasibility of 
macroscopic quantum tunneling in the intrinsic Josephson junctions of a high temperature super- 
conductor Bi2Sr2CaCu208+a, and find it to be characterized by a high classic-to-quantum crossover 
temperature and a relatively weak quasiparticle dissipation. 
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A marked feature characterizing cuprate high temper- 
ature superconductors (HTSC) is its strong two dimen- 
sionahty. In particular the bismuth-based HTSCs, for 
which this anisotropy is prominent, are best viewed as 
stacks of superconducting Cu02 planes weakly linked 
through intrinsic Josephson junction (IJJ) type couplings 
0. These built-in atomic scale links are taylor-made for 
technical applications difhcult to achieve with artificial 
Josephson junctions (JJ), and its study has now devel- 
oped into an active interdisciplinary field. We report 
below what is to our knowledge the first successful ob- 
servation of the macroscopic quantum tunnelling (MQT) 
0, 0, 0, 13 of the phase variable of the superconducting 
order parameter through the potential barrier of an IJJ, 
opening up an entirely new direction for HTSC appli- 
cations. While the corresponding phenomena had been 
observed at around 300 mK in conventional J Js |^ |H , 
we have confirmed MQT behavior at approximately 1 K 
apparently reflecting the characteristically high plasma 
frequency of IJJs. Our results are highly nontrivial in 
that they also demonstrate the feasibility of MQT in spite 
of the presence of dissipative low energy quasiparticles 
mm [13, which is the other hallmark of HTSCs. 

Current-biased JJs offer an ideal stage for realizing a 
variety of macroscopic quantum phenomena, e^. energy 
level quantization within the potential well |j,llllll3l and 
the associated MQT and macroscopic quantum coherence 
0,0, all of which have come to be recognized as having 
immediate implications for qubit applications Emilia 
Ir^ . In particular, a phase qubit utihzing MQT has been 
reported by Martinis and co-workers d5]- 

Aside from external noises and disorder, a primary 
source which stands as an obstacle towards observation 
of MQT is the influence of non-superconducting quasi- 
particle excitations In conventional s-wave super- 
conductors, all quasi-particle states are separated from 
the superconducting ground state by a finite energy gap 
and thus become essentially inaccessible upon going to 
sufficiently low temperatures; hence the observability of 



MQT. The situation is drastically altered when wc turn 
to HTSCs, which are c?-wave superconductors. The latter 
are characterized by four nodes in the order parameter 
at which the energy gap vanishes [ilEl. HTSCs there- 
fore necessarily sustain a finite population of dissipative 
quasiparticles down to the lowest temperatures, which 
would lead one to render MQT unattainable. While this 
accounts for the absence of previous phase-MQT reports 
in HTSCs, several theoretical works have recently put 
this naive view into question by inferring that dissipation 
due to nodal quasiparticles is not strong enough to totally 
destruct the coherence required for MQT l7i|^|9|,il0j. In 
addition to the moderately weak influence of quasiparti- 
cles, it is also crucial from the experimental point of view 
to have JJs bearing a high degree of underdamping 
and hysterics, and ||2Jl a measurement system sufficiently 
detached from the environment. Devising samples and 
apparatus fulfilling this prerequisite constituted an inte- 
gral portion of our project of observing MQT in the IJJs 
of BiaSraCaCuaOg+a (Bi-2212). 

We now turn to the description of our samples and 
measurement set-up. In Bi-2212, it is well known that 
the stacking crystal structure along the c axis consists 
of an alternating array of the superconducting (S) CuOa 
layers and the insulating (I) BiO and SrO layers. The 
excellent JJ properties exhibited by IJJs in Bi-2212 is 
by now well established; being ideal atomic sized junc- 
tions P, [23, they are in particular free of defects and 
surface roughness of the insulating oxide layers which 
degrade the quality of artificially made JJs, such as in 
Al and Nb. To probe the characteristics of the built-in 
IJJs, it is necessary to construct a pass for the interlayer 
current flow as illustrated in Fig. 1(a) employing the 
combination of the double side etching and 3-D focused 
ion-beam etching techniques 0,113 . The junction thus 
fabricated will generally consist of several layers of IJJs 
coupled inductively and via charging effects -the num- 
ber of which is technically difficult to control owing to its 
minute (Angstrom) scale. Therefore, the IJJs form a ID- 
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barrier whose rate is given by [26 
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FIG. 1: (a) Scanning ion-beam micrograph of the IJJs. The 
critical temperature of the sample is around 90 K. The num- 
ber of elemental junctions in the IJJs is less than 20 and the 
junction dimension in the ab plane is 1.76 x 0.86 fim^ . (b) 
Schematic of the bias lead configuration. Electrical properties 
of IJJs were measured with the four probes method, (c) I-V 
characteristic of the IJJs measured at 4.2 K. The I-V curves 
were measured by ramping a bias current up and down re- 
peatedly. A total of two resistive branches (first and second 
branches) are exhibited in this figure. The inset shows V or 
/ vs. time. 



array of S-I-S type J J. Close to zero bias, the behavior is 
nevertheless well captured in terms of a single junction 
model |23l |. as evidenced from our experiments detailed 
below. The IJJs were mounted on the copper-mixing 
chamber of a '^He/'^He Oxford dilution refrigerator with 
an excellent thermal contact and cooled down from 4.2 K 
and 40 mK. All cable connections at room temperature 
were guided through high impedance cupronickel semi- 
rigid coaxial cables, which work as a low-pass filter above 
1 GHz, to avoid propagating thermal flows and external 
noises into the sample. Moreover, by using LC low-pass 
filters along these cables we achieved a strong attenuation 
of the external noises above 15 MHz. The dilution refrig- 
erator along with the entire analog-measurement system 
was placed in a shielding room. 

The dynamics of a Joscphson junction [3 is con- 
veniently described as the motion of a fictitious parti- 
cle in a washboard potential C/(0) = —Uo{x(p + cos(f>), 
where the phase difference </> is to be viewed as the par- 
ticle coordinate. Here Uq = hlc/'2e is the Josephson 
coupling energy and x = I/Ic {I- bias current, Ic- criti- 
cal current). The particle at the bottom of the poten- 
tial well oscillates at the plasma frequency defined as 
u>p = ujpo{l — x^)^!^, where Wpo — {2elc /hCY/'^ is the 
zero-bias plasma frequency and C is the junction capac- 
itance. In the presence of a bias-current, U{<t)) consists 
of a periodic potential structure separated by the energy 
barrier £;(x) = C7o{-7rj:-t-2[xsin"^ a;+(l-a;2)]i/2}. Ther- 
mal fluctuations at sufficiently high temperatures assists 
a thermally-activated escape of the particle across the 
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where fcs is the Boltzmann constant and T is the tem- 
perature. As T is lowered, thermal activation is sup- 
pressed exponentially and the escape process eventually 
gets dominated by MQT events. This crossover takes 
place at the temperature T* ~ hLUp/2TrkB 0,113, and 
the MQT rate is expressed as 
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The current-voltage [I-V] characteristic of the IJJs 
measured by biasing along the c axis (see Fig. 1(b)) 
at 4.2 K is displayed in Fig. 1(c). The typical I-V char- 
acteristic of IJJs exhibiting a multi-branch structure was 
observed. The total number of branches corresponds to 
the number of elemental junctions in the IJJs as one can 
confirm by applying a larger bias. When biased, a switch- 
ing to a finite voltage state takes place at a random value 
of the switching current Isw and the I-V curve shows a 
hysteresis. The zero-bias Josephson current is slightly 
larger than the first branch of IJJs, which indicates that 
the noise level of our system is quite low. In addition, the 
McCumber parameter (3c=2eIcC / h and the damping 
parameter Q~^={ujpoCR)~^ were ~ 10^ and ~ 10~^, re- 
spectively. R is the effective shunt resistance, and Rq is 
the quasiparticle resistance which is estimated from the 
slope of the first branch of an I-V curve at zero bias. 
This high degree of hysteresis and the underdamping are 
crucial for detecting the MQT. 

The physical nature of the escape process in JJs can 
be elucidated by analyzing the statistics of the stochas- 
tic switching from a zero-voltage state to a 
finite voltage state. To this end measurement of Isw 
was repeated for 2000 times at each temperature with a 
current resolution A/ of 20 nA, which we checked were 
sufficient to determine the switching current distribu- 
tion P{Isw)- Figure 2(a) shows the temperature depen- 
dence of P{Isw)- The measurements were performed 
for the first branch of IJJs. In this experiment a ramp- 
shaped waveform was used with the constant bias speed 
dl/dt. Histograms of P{Isw) were calculated against 
the bias current and the bin width corresponds to the 
current resolution A/ of the measurement setup. As T is 
lowered, P{Isw) narrows and shifts to higher currents. 
The P{Isw) measured at 4.2 K are fitted with the ther- 
mal activation theory proposed by Kramers 0, |0 , and 
the zero-noise critical current obtained from the fitting 
is /c=48.54±0.02 /iA. This value gives a practical esti- 
mate of 7(7 at K. The results exhibit good agreements 
with the Kramers theory at 4.2 ~ 1.1 K, but deviates 
substantially below IK. Note also that below IK, the 
amount of shift of P{Isw) along the current axis reduces 
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considerably. These findings strongly infer that quan- 
tum tunneling takes over and the escape process is es- 
sentially T-independent in this regime as verified below. 
We also expect that the good fitting of P{Isw) in the 
thermally activated regime to conventional theories is an 
indicator that the weakest link in our IJJs effectively be- 
haves like a single Josephson junction. Figure 2(b) shows 
T~^(Iqur) calculated from the experimental data in Fig. 
2(a) miiil. As T is lowered, the slope of In t-^{Isw) 
increases and t~^{Isw) shifts to higher Isw similar to 
the behavior of P(Isw)- One observes from a compari- 
son of the two figures that below 1 K, where the amount 
of shift of P{Isw) gets smaller and saturates, the rates 
T~^{Isw) overlap with each other. This implies that 
T~^{Isw) is also independent of T in this regime. Using 
the value of Ic obtained from the fitting of P{Isw), the 
tunneling rates are estimated from Eqs. iQJ and 10). The 
fitting of T^^{Isw) agrees with the experimental data. 

The standard deviation a(T) of P{Isw) is plotted 
against T on a double logarithmic scale for two sam- 
ples in Fig. 3. In the high T region (1-4.2 K), cr(T) 
is proportional to T^/^ as expected, and decreases with 
decreasing T owing to the suppression of thermal activa- 
tion. But below 1 K it is saturated and independent of 
T in accordance with Fig. 2. T^gc |^ also exhibits a T- 
independent behavior, (see the inset of Fig. 3) Actually, 
we have also observed this kind of saturating behavior of 
cr(T) in the IJJs fabricated in Bi-2212 whisker single crys- 
tals. These phenomena agree with the tendencies which 
have been reported for low- Tc J Js 0, . To check that 
the observed saturation is not due to external noise, we 
confirmed that connecting LC low-pass filters in series 
does not effect the value of a{T). Moreover, an inter- 
mittent ramp bias was applied to the sample in order to 
avoid heating. We are thus lead to conclude that the ex- 
perimental results presented here are direct consequences 
of MQT of the macroscopic junction phase of IJJs fab- 
ricated in HTSC Bi-2212. The effect of MQT appears 
below the crossover temperature T*, which we estimate 
from the data to be 0.75 K for Jc=48.54 ^A and 0.95 K 
for /c=84.92 /zA, which are shown in Fig. 3 by the open 
and filled arrows, respectively. Meanwhile, based on the 
uncertainties of Ic and C, and the bias dependence of T* 
as obtained from the formula T* — hujp/2TrkB 0,113 is 
0.75~0.80 K for /c=48.54 and 0.85-1.0 K for 84.92 
fiA. These results agree well with the experimental values 
of T*s. 

Having established the occurrence of MQT, let us 
briefly account on its physical implications. We begin 
by mentioning two factors essential to our findings, i.e. 

the large magnitude ujp of the overall energy scale 
involved, and (01 the relatively weak form of dissipation 
with which the quasiparticle transport affects the system. 
The former aspect is apparently reflected in the thermal- 
to-quantum crossover temperature, which to the best of 
our knowledge is the highest reported to date. Though 
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FIG. 2: (a) Switching current distribution P(Isw) of IJJs. 
The solid and dashed lines are the fitting |2a.l29l| in the ther- 
mal and MQT regions, respectively, and the zero-noise critical 
current Ic was extracted from these fittings. The measure- 
ment condition and the junction parameters were the follow- 
ing: dI/dt=42A mA/sec, A7 = 20 nA, /c=48.54±0.02 /xA, 
C=76.26 fF, ujpo/2-K=221.3 GHz. The junction capacitance 
C was estimated by C = er£oS/d, where d = 15 A P, l20l| is 
the c-axis lattice constant of Bi-2212, S is the junction area, 
Er = 10 20] and £o are the relative permittivity and the vac- 
uum permittivity, respectively, (b) Escaping rate t~^{Isw) 
of IJJs. The theoretical thermal rates calculated by Eq. (0 
are shown as solid lines, and MQT rate without dissipation, 
which was calculated by Eq. (|5J is exhibited as dashed line. 

the temperature itself is not the main emphasis here, it 
does expose a potential advantage enjoyed by HTSC ma- 
terials, which deserves to be explored further. Meanwhile 
the latter factor concerns the much more subtle issue 
of "taming" the nodal quasiparticles mentioned earlier, 
which bring us back to the I-V characteristics displayed 
in Fig. 1(c). In general, the quasiparticle dynamics in 
a JJ is manifested in the return curve (the low-voltage 
subgap portion of the curve following the jump to the 
first branch, obtained upon switching off the bias cur- 
rent), a regime where Cooper pairs cannot participate 
in the charge transport. In J J stacks artificially fabri- 
cated using conventional superconductors, this curve ex- 
hibits an abrupt drop to nearly zero current at low tem- 
peratures |20l |. signalling a strong suppression of quasi- 
particle excitations. Meanwhile, the resistively-shunted- 
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cross- fertilization between the fields of HTSC science and 
quantum information. 
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Note added- After completion of this work, we 
learned that Bauch et al. have also observed the MQT at 
around 40 mK in an HTSC YBa2Cu307_5 grain bound- 
ary Josephson junction js^- 



FIG. 3: Standard deviation a{T) of P{Isw)- <y{T) were cal- 
culated by £7 = {{I%w) ~ {IswYY^"^ for two samples which 
had different Ics- The solid lines show the theoretical fit- 
ting of the thermal region. a{T) starts to saturate below 
the crossover temperature T* [12 . l27l . The experimental T*s 
of the measured samples are around 1 K, falling within the 
shaded crossover region calculated theoretically. The inset 
shows the escape temperature Tcsc vs. T. The dotted line 
corresponds to T=Tesc. 

junction model widely employed in quantum mechanical 
treatments of dissipation in JJs 0] can lead to the domi- 
nance of ohmic behaviour and decoherence when the re- 
sistance becomes a relevant perturbation. The power- 
law-like tail seen in our data falls into a novel class of 
behaviour termed the superohmic dissipation |3Cll |. which 
is intermediate between these two cases. That this sub- 
tle form of dissipation should play a crucial role in our 
observations is in accord with some of the recent theories 
on d-wave J Js HH H [l^ ■ 

Our study bores out a number of fundamental ques- 
tions, which reside on the intersection of two major areas 
of current research: HTSC physics and future quantum 
technology. For example it provides a fresh perspective 
on the issue of differentiating between very distinct forms 
(e.g. transverse-momentum conserving and nonconserv- 
ing processes) of quasiparticle tunnelling between d-wave 
superconductors, several of which may account for super- 
ohmic dissipation. Another related problem -with direct 
implications to optimizing device geometry- is to resolve 
how the differences between generic d-wave JJs and the 
interlayer IJJs would influence MQT experiments. For 
instance, the absence of zero-energy bound states H^and 
the presence of a bilayer tunnelling matrix element sup- 
pressing nodal quasiparticle tunnelling jsj] in the latter 
can both be considered advantages peculiar to the IJJs 

In summary we have successfully observed phase MQT 
in the IJJs of HTSCs, inviting further investigation 
on coherence-dissipation competition in these systems. 
Though yet an initial step towards applications, we an- 
ticipate that efforts along this line may lead to a rich 
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